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1. Introduction 



The Hubbard model has been studied as the prime example of a microscopic theory of 
strongly correlated electrons for several decades. More recently it was proposed as a model for 
superconductivity by P.W. Anderson. It seems to be clear that in order to make close contact 
to experimental results, the Hubbard hamiltonian has to be modified by adding competing 
interactions. Despite the simple form of its hamiltonian, the mathematical structure of 
the Hubbard model is rather complex (as can be seen in the Bethe Ansatz solution in 
one spatial dimension) and exact results about its physical features are quite difficult to 
obtain. The ground states in more than one dimension are still unknown, and the question 
under what conditions the model exhibits superconductivity is yet unanswered. Therefore 
it would be useful to have a model similar to the Hubbard model, where the physics is more 
transparent. In a recent publication we have presented a new model, describing strongly 
correlated electrons on a general d-dimensional lattice'^', which differs from the Hubbard 
model by moderate interactions of nearest neighbours. The nature of these interactions is 
very similar to the ones of Hirsch's model of hole superconductivity '^l In one dimension we 
the model is integrable. We gave an exact expression for the ground state in the attractive 
regime in 1, 2 and 3 dimensions and showed that the ground states of the model for attractive 
and moderately repulsive on-site interaction exhibit off-diagonal long range order (ODLRO) 
and are thus superconductingt^l . In order to gain further insights into the physical features 
of the model it is reasonable to first study its one dimensional version, where integrability 
permits us to obtain exact results. In this paper we construct the nested Algebraic Bethe 
Ansatz (NABA) for the model and use it to determine higher conservation laws. We also 
prove a lowest weight theorem for the Bethe Ansatz states. 

In two forthcoming publications we will determine the structure of the low lying excitations 
over the ground states and prove completeness of the ^(21 2) extended Bethe Ansatz W. 

Electrons on a lattice are described by operators ^ , j = 1, . . . , L, cr = ±1, where L is 
the total number of lattice sites. These are canonical Fermi operators with anti-commutation 
relations given by {ct^., c^.,.} — 5j^ j5^.^. The state |0) (the Fock vacuum) satisfies c^^^\^) — 0. 
There are four possible electronic states at a given lattice site i 



|o), IT)i = 4,|o), U), = 4_,|o), |iT), = 4_,4,|0) . (1.1) 

By Hj = A,(T%a denote the number operator for electrons with spin a on site i and we 
write n^^n^ ^+n^ _^. 

The hamiltonian on a lattice of L sites in the grand canonical ensemble is 

L L L 

H^H° + Uj2inj,i - 2)K-i - 2) - ^E^.- + ^^EKi - ^.,-1) (1-2) 
j=i j=i j=i 
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where = - Y^^j=i ^^j+i given by 

+ 2 (^i - - 1) + + S-,-iSm4+i,ic3+i,-i 

11 11 

+ Ki - 2)K-i - 2^ + ~ a^^^'^'+i'-i ~ 2^ • 

(1.3) 

Here U is the Hubbard model couphng constant, h is an external magnetic field, and ji is 
the chemical potential. 

It was shown in [1] that is invariant under a u(2|2) symmetry algebra. The three other 
terms in H are elements of the Cartan subalgebra of m(2|2) (and commute with and 
each other), and thus and H have a complete set of simultaneous eigenstates. In what 
follows we will construct a complete set of eigenstates of (and of H) by using a nested 
Algebraic Bethe Ansatz and the m(2|2) symmetry of H^. The algebra u{2\2) has a total of 16 
generators, 8 of which are fermionic and bosonic, respectively. The eight bosonic operators 
fall into two su{2) and two u{l) subalgebras. The two su{2) algebras are generated by the 
spin-operators [We shall always give local expressions O • for symmetry generators, implying 
that the global ones are obtained as C = ^^=1 Oy] 

and by the ry-pairing like operators'^' 

Note that unlike for the case of the Hubbard model as treated in [5] there is no factor of 
(— 1)^ in the definition of rjj. The two u{l) charges are given by the identity operator and 
the operator 

^ = E^.- ^, = Ki-^)K_i-^) . (1.6) 

As all eigenvalues of + 7]^ — —N_i + ■§ (where 7V_^ = number operator 

for spin down electrons) are either integer (for even L) or half-odd integer (for odd L), 
the two 5^(2) algebras are not quite independent. Although all six generators of the two 
su{2) algebras commute with H^, the Hilbert space of eigenstates of will only carry 
representations of su{2) x su{2) with integer (for even L) or half-odd integer (for odd L) 
quantum numbers. The same kind of restriction carries through to the complete ^(212) 
symmetry. The eight fermionic generators are given by 
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and their hermitean conjugates. 

We have shown in [1] that can be written as minus the sum over graded permutation 
operators 

L 

//° = -^n^'^+^ , (1.8) 

where the operator n-''^+^ permutes the four possible configurations (1.1) between the sites 
j and j + 1, picking up a minus sign if both of the permuted configurations are fermionic, 

n^'^-+i|o),. X |o),.+i = |o),. X |o),.+, 
n^'^'+^UT),xUT),+i = |iT),x|iT),+i 

n^'^+i|0),. X = \a)^ X |0),.+i (1.9) 

etc. 

It is clear that this form of interaction conserves the individual numbers A^^^ and A^j^ of 
electrons with spin up and spin down, and the numbers A^; and A^;^ of doubly occupied 
("local electron pairs") and empty sites ("holes"). We will choose the following conventions 
throughout this paper 

= number of single electrons with spin up 

— number of single electrons with spin down 

— + — number of single electrons 
Ni — number of local electron pairs 

Nf^ — number of holes 

Nf, — Nfj^ + Ni — number of "bosons" . 

The outline of this paper is as follows : 
In section 2 we perform a detailed construction of the Algebraic Bethe Ansatz of the model. 
We derive six different forms for the Bethe Ansatz Equations (BAE) and the eigenvalues 
of the transfer matrix. The graded Quantum Inverse Scattering Method (QISM), discussed 
in section 2, enables us to obtain expressions for the hamiltonian and (an infinite number 
of) higher conservation laws at the quantum level. These conserved charges are of interest, 
because physical interactions, although of short range, are not generally well approximated 
by interactions involving only nearest neighbours. The charges under consideration involve 
interactions of longer range (next nearest neighbours, next next nearest etc.) and can be 
added to the hamiltonian without spoiling the integrability of the model. Thus it is possible 
to construct integrable models with longer range interactions by using higher conservation 
laws'^1 . The first nontrivial higher integral of motion is for example given by the expression 

L 

i^(3) = ^5][i^,Vl>^^l,.] > (1-10) 

k=l 
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where -f^^ is the density of the hamiltonian defined in (1.2). Section 3 is devoted to the 
derivation of exphcit formulas for higher conservation laws. 

2. Graded Quantum Inverse Scattering Method 

In this section we discuss how to embed the hamiltonian (1.3) into the framework of the 
Quantum Inverse Scattering Method (QISM). For supersymmetric theories it is necessary 
to modify the QISM along the lines discussed by Kulish and Sklyanin in [7,8]. In [9] we 
gave a summary of the "graded" version of the QISM, so that we will constrain ourselves 
to a very brief discussion below and refer to [9] for more details. Starting point of the 
QISM is an it!-matrix, which obeys a graded Yang Baxter equation. Prom this i?-matrix 
we then construct a family of commuting transfer matrices, describing a "fundamental" spin 
model. The transfer matrix is the generating functional of an infinite number of conservation 
laws, the second of which we show to be the hamiltonian (1.3). Finally we construct a set of 
simultaneous eigenstates of the transfer matrix and the hamiltonian, using a nested Algebraic 
Bethc Ansatz (NABA)'^'-'^'^"^^'^! . Due to the grading there exist a priori 4! = 24 choices for 
the i?-matrix, all of them describing the same physical system, but leading to different (yet 
equivalent as shown in Appendix B) forms of the NABA. As the model has the two trivial 
discrete symmetries of "spin refiection" (interchange of spin up and spin down electrons) 
and of interchange of empty sites and doubly occupied sites, there are essentially only six 
different solutions^. One of the six solutions has been previously obtained in [7]. By the 
same kind of reasoning one can show that there are ^"j^"^)' Bethe Ansatz solutions for a 

^ n! ml 

permutation-type model with M(m|n) symmetry. 

2.1. Yang Baxter Equation 

Consider the graded linear space V^"^^^\ where m and n denote the dimensions of the 
"even" and "odd" parts respectively. A basis {e-,^, . . . , e^^,„} of \/(™'+") consists of m + n 
vectors with Grassmann parities such that for m vectors = and for n vectors 
= 1. The grading of the matrix elements M^^ of linear operators M on (in the 

basis {e;^, . . . , e^_,_„}) is given hy ej^ ^ = e^-\- e^,. The supertrace of M is defined as 



The graded tensor product space (g) is given in terms of its basis vectors 

{^a ® b = 1-1 . . . iTn -\- n} &s follows 




(2.1) 



a=l 



e„®efe)v„ti;b(-l)'''a 



(2.2) 



^Each of these solutions generates a "multiplet" of four NABA via the substitutions N^^ A''j^ and Ni <-> Nj 



h 
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The action of the hnear operator F (8) G on the vector v<Siw in (g) is defined by 

{F^G){v®w) ^ F{v)^G{w) , (2.3) 

which results in matrix elements of the form 

(^®G'fi = F„,G',,(-l)^c(^<'+^^) . (2.4) 

The identity operator / in I/*^™'"-* and the matrix 11 that permutes the individual 

linear spaces in the tensor product space are given by 



U(v®w)^(w®v), (n)2^=5„^,/,^,^(-i)^H% . 

A matrix -R(A) (depending on a spectral parameter A) is said to fulfill a graded Yang- 
Baxter-equation, if the following identity on F^"*!") O V^""^""^ (g) F^"*!") holds 

(1 R{X - fi)) {R{X) ^ I) {I ^ R{li)) ^ {R{l^) ^ I) {I ^ R{X)) {R{X - 1^) I) . (2.6) 

In components this identity reads 

Note that despite the fact that the tensor product in (2.6) carries a grading, there are no 
extra signs in (2.7) compared to the nongraded case. It is easily checked that the R-matrix 



a(A) = -A_, 6(A) ' ^^-^^ 



R{X) =6(A)/ + a(A)n 

X+i ^ ' x+% 

fulfills equation (2.7). 

2.2. Construction of the Transfer Matrix 

Prom (2.7) one can derive the equation 

R,,{X - fx) m,,R,,{X)] ® [U,,R,M]) = {[U,,RM] ® [U,MX)]) R,,{X - /x) , (2.9) 

where the indices 1,2,3 indicate in which of the spaces in the tensor product space 

y(m|n) ^yiiiiin) ^ y(m|n) ^j^g matrices act nontrivially. The tensor product in (2.9) is between 
the spaces 1 and 2. We now call the third space "quantum space" and the first two spaces 
"matrix spaces" . The physical interpretation of the quantum space is as the Hilbert space 
over a single site of a one-dimensional lattice. We now consider the situation, where inter- 
twining relations of the type (2.9) hold for all sites of a lattice of length L. The quantum 
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space index "3" now gets replaced by an index labelling the number of the site. We define 
the L-operator (on site n) as a linear operator on KnaJril (where TC^ — l/^"*'") is the 
Hilbert space over the n*'* site, and V^^}^ is a matrix space) 

LnW% = ^ZR^>^)1^ = (KA)n + a{X)ir^^ . (2.10) 

is a quantum operator valued (m + n) x (m + n) matrix, with quantum operators acting 
nontrivially in the n^^ quantum space (of the direct product Hilbert space over the complete 
lattice (8)j=i'Wj)- The greek indices are the "quantum indices" and the roman indices are 
the "matrix indices". Equation (2.9) for the n*'* quantum space can now be rewritten as the 
operator equation 

R{X - I,) (L„(A) ® LM) = (LM ® ^n(A)) R{X - ■ (2.11) 

Here the graded tensor product is between the two matrix spaces and R only acts in the 
matrix spaces. The intertwining relation (2.11) enables us to construct an integrable spin 
model as follows. 

We first define the monodromy matrix Tj^{X) as the matrix product over the L-operators on 
all sites of the lattice, i.e. 

TM = LML,_,{\)...L,{\) 
((r,(A))-) . . = L,{\r:^^,L,^MCj,_, ■ ■ ■ L,iX)t,, X 

P1---PL 

X ( — l)S3=2(^«j+'^/'3)Si=l^"i 

T^(A) is a quantum operator valued (m + n) x (m + n) matrix that acts nontrivially in the 
graded tensor product of all quantum spaces of the lattice and by construction fuUfiUs the 
same intertwining relation as the L-operators (as can be proven by induction over the length 
of the lattice) 

R{\ - ii) (TJA) ® TM) = {TM ® T^W) R{X - t^) ■ (2-13) 

The transfer matrix t{X) of the integrable spin model is now given as the matrix supertrace 
of the mondromy matrix 

m+n 

r(A) = str{T,{X)) = j;(-l)^''(T^(A))- . (2.14) 

a=l 

As a consequence of (2.13) transfer matrices with different spectral parameters commute. 
This condition implies that the transfer matrix is the generating functional of the hamiltonian 
and of an infinite number of "higher" conservation laws of the model. 
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(2.12) 



2.3. Trace Identities 



Taking logarithmic derivatives of the transfer matrix at a special value of the spectral 
parameter, one can generate higher conservation laws'^^^ . For the transfer matrix constructed 
from the i?-matrix (2.8) along the hnes discussed in section 2.2., the corresponding hamilto- 
nian is obtained by taking the first logarithmic derivative at zero spectral parameter 



. d\0g{T{X)) 



= -^(n'^'^^^^ - 1). (2.15) 

A=0 k=l 



The proof of this identity can be carried out in the same way as for the ungraded case, 
the main difference being the grading of the tensor product of the quantum spaces (see 
(2.12)). By shifting the energy eigenvalues by a constant we obtain the expression (1.8) for 
the hamiltonian (1.3) 



^o^_.^log(r(A)) 



dX 



A=0 



L ^ H^,^-L , (2.16) 



if we choose our underlying graded vector space to have signature (2, 2), i.e. to have a basis 
with two fermionic and two bosonic states. This shows that the transfer matrix constructed 
from the L-operator (2.10) and i?-matrix (2.8) is indeed the correct transfer matrix for the 
hamiltonian (1.3). Higher conservation laws are obtained as the coefficients of the power 
series 

log(r(A)(r(0))-i) = . (2.17) 

fc=i 

There exists however a simpler method for the construction of higher integrals of motion 
than taking logarithmic derivatives, which we will discuss in section 3. For fundamental 
spin models (and thus for the model at hand) the transfer matrix at the special value of the 
spectral parameter (which is zero for our case) is equal to the translation operator. Therefore 
we can construct the momentum operator from the logarithm of the transfer matrix at zero 
spectral parameter 

P = log(T(0)) . (2.18) 

By construction the momentum operator commutes with the hamiltonian and all higher 
conservation laws i.e. [P, — — [P, ifj-fc+i)]. 

2.4. Algebraic Bethe Ansatz with a bosonic background (FFBB grading) 

Due to the fact that there are four different configurations per site for the model defined 
by (1.2), the Hilbert space at the k^^ site of the lattice is isomorphic to (D^ and is spanned 
by the four vectors = ( 1 0)^, 62 = (0 1 0)^, 63 = (0 1 0)^, and 
64 = ( 1 )'^. In this section we consider a grading such that and 63 are fermionic 
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(representing spin down/spin up electrons respectively) and and are bosonic (doubly 
occupied/empty site). In terms of the Grassmann parities this means that = €2 — 1 and 
^3 — Q — 0. We pick the reference state in the A;*'* quantum space |0);;., and the vacuum |0) 
of the complete lattice of L sites, to be the bosonic completely unoccupied state, i.e. 



|0)n 








|0) = ®ti|0). 



(2.19) 



This choice of grading implies that the L-operator (2.10) is given by the following expression 



/a{X)-b{X)ei' 




b{\)ef 


b{\)ef 




-KA)e^^ 


a(A) - b{X)ef 


b{X)ef 


b{X)ef 






h{X)ef 


a{X) + b{X)ef 


b{\)ef 




V KA)ei^ 




b{X)ef 


a{X) + b{X)e'^ 


I 



(2.20) 



where e^^ are quantum operators in the n^^ quantum space with matrix representation 
{^'n)ai3 — ^aa^bi3- The monodromy matrix (2.12) is a quantum operator valued 4x4 matrix, 
which we represent as 



T^{X) = L^{X)L^_,{X)...L,{X) 



fA,,{X) A,(A) A,,{X) B,{X)\ 

A,,{X) A,,{X) A,,{X) B,{X) 

A,,{X) A,,{X) A,,{X) 53(A) 

\C,{X) C,{X) C,{X) D{X)/ 



(2.21) 



The transfer matrix is then given by 

r{fj,) = str{TM) = -AM - A2(/^) + AM + D{y) 



(2.22) 



The action of Lj^[X) on the reference state on the A;*^ site can be easily determined by using 
the matrix representation of the e'^ in (2.20) and the vectorial representation of \Q)f^ (2.19) 



( o(A) 0^y 

a(A) 

a(A) 

V6(A)e^^ b{X)ef b{X)ef \) 



|0). 



(2.23) 



Using (2.12) and (2.23) we determine the action of the monodromy matrix on the reference 

state to be , , . . r 

({a{X)f 0\ 

(a(A))^ 

(a(A))^ 

V Ci(A) C^iX) C3(A) 1/ 



7^l(A)|0) = 



(2.24) 
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We will now construct a set of eigenstates of the transfer matrix using the technique of 
the NABA. Inspection of (2.24) reveals that Cj(A) are creation operators with respect to our 
choice of reference state. This observation leads us to the following Ansatz for the eigenstates 
of t(/x) 

|A„...,AJF) = C,^(AJC,^(A2)...C„JAJ lOF'^--^ , (2.25) 

where the indices run over the values 1,2,3, and F"" ' '*i is a function of the spectral 
parameters A^-. The action of the transfer matrix on states of the form (2.25) is determined 
by (2.24) and the intertwining relations (2.13). The components of the intertwining relations 
(2.13) needed for the construction of the NABA are 

-(-^'-'"^^'("'^"W (2.20) 



where 



CaMi) C„,(A2) = r{X, - X,)li:i C,^(X,) C,^(X,) 



(2.27) 



Here H^^^"^ is the 9x9 permutation matrix corresponding to the grading = eg = 1 , = 0. 
r{fj,) can be seen to fulfill a (graded) Yang-Baxter equation on its own 

- H)ZZ rintil = rinralU r{X):itl r(A - , (2.28) 

and can be identified with the i?-matrix of a fundamental spin model describing two species 
of fermions and one species of bosons. Using (2.26) we find that the diagonal elements of 
the monodromy matrix act on the states (2.25) as follows 



D(/.)|A„ . . . , AJF) = n ^(^^ _^) |Ai, . . . , A„|F) 



J 

n 



k=i "1-"" j=i 

{-AM - AM + ^33(/"))|Ai, . . . , AJF) = 

n ^ n 



(2.29) 



(2.30) 



.ai 



k=l 3=1 
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where 



■K 

..a„ 



str(LW(/x-AjL«,(/x-A„_J 



Ai)) 



and 



6(A)n(i) + a(A)/(^) 
/a(A) 

n(i) r(A) = 



(2.31) 
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^31 



v 



-b{X)ef b{\)e: 

-b{X)el? a(A) - bi\)ef b{\)ef 

b(X>^? b{X)e'^ a(A) + 6(A)e: 

Note that the tensor products between the n quantum spaces in the expression for the 
operator t^^\ii) are once again graded, which resuhs in minus signs given by (2.12). L« 
and r{n) can be interpreted as L-operator and /^-matrix of a fundamental spin model (r 
fulfills the Yang-Baxter equation (2.28)), describing two species of fermions and one species 
of bosons. Hence Tn^ (/i) and t^^^ (/x) are the monodromy matrix and transfer matrix of the 
corresponding inhomogeneous model. This model can be identified with the inhomogeneous 
supersymmetric t-J model on a lattice with n sites. The Algebraic Bethe Ansatze for the 
supersymmetric t-J model were constructed in detail in [9], so that we will restrict ourselves 
to a abbreviated discussion and refer to [9] for derivations of some of the results. Inspection 
of (2.29) and (2.30) together with (2.22) shows that the eigenvalue condition 



r(/x)|\,...,AjF) = K/^) |Ai,...,A„|F) 



(2.32) 



leads to the requirements that F ought to be an eigenvector of the "nested" transfer matrix 
T^^\n), and that the "unwanted terms" cancel, i.e. 



(2.33) 



The notation for the eigenvalues of the transfer matrix r(//) is to be considered as a 
shorthand as the eigenvalues depend also on the spectral parameters A^ and coefficients F. 
The quantities h.^. and A^^. are computed in Appendix A. Using their explicit expressions in 
(2.33) we obtain the following conditions on the spectral parameters A^ and coefficients F, 
which are necessary for (2.32) to hold 



n 



(2.34) 



This completes the first step of the NABA. In the next step we will now solve the first nesting. 
The condition that F ought to be an eigenvector of t^^\ii) requires the diagonalisation of 
r''^'*(^), which can be carried out by a second, "nested" Bethe Ansatz. From (2.28) and 
(2.31) the following intertwining relation is easily derived 



r(A - /.) (tW(A) rW(/.)) = {T^'\^Ji) T(^)(A)) r(A - 



(2.35) 
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If we write 

Ti'\f^) = 41V) 42V) , r('\f^) = - 42V) + ^^^H/") , 

(2.36) 

then (2.35) and (2.27) imply that 

^^i"""^^''^' (2.37) 

Here the roman indices take the values 1,2, both of which are fermionic (e^ = 1 = 63). The 
i?-matrix r^^\fj,) can again be proven to fulfill a Yang-Baxter equation of the form (2.28). 
As the reference state for the first nesting we pick 

|0)i^)= j , |0)«=®Li|0)i'^ . (2.38) 

The action of the nested monodromy matrix Tn^^ {/j.) on the reference state |0)^^^ is determined 
by (2.31) and we find 



(2.39) 



D«(^)|0)« = |0)« . 

We now make the following Ansatz for the eigenstates of r*^^^ (//) 

|A«,...,A«|G) = cW(aW) (A(^))...cW(A«) |0)« GS-^. . (2.40) 

Here the indides can take the two values 1 and 2. These states can be related to the 
coefficients in the following way : The state jA^^^ . . . , Xn^) "lives" on a lattice of n 

sites and is thus an element of a direct product over n Hilbert spaces. In components it 
reads \X\^^ . . . Xn})„ „ , which can be directly identified with 

The action of r*^^^(//) on the states (2.40) can be evaluated with the help of (2.37) and (2.39) 



"^1 , 

D«(/.)|Ai^ . . . , AW|G) =n — ti^-IaS^ . . . , A«|G) 



J 

n 



k k "• (2-41) 

+ E (Ai") '■' n C«(a5'>)|0><')G"". «. . 



fe=l "1- -"1 j=l 
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(-A^S(^^)-A^^(^^))\x?,...,xS\G)^ 

1=1 j=i ^i/^' - ) 1=1 (2.42) 



+ 



fe=l ^ "1 7 = 1 



..ai 



Here r^^^ {ji) is the transfer matrix of the second nesting 

= ML^li^. - A«)Lif_,(;. - A«_,) . . .LrV - A^^^)Lf)(^ - A«)) , 

where 

(A) = 6(A)n(^) + a(A)/(2) 

V -KA)e^^ a{\)-b{X)ef) 

(2.43) 

The operator T'^^^(/Lt) can be interpreted as the transfer matrix of an inhomogeneous model 
on n-^ sites describing two species of fermions. It is easily seen from (2.42) that for the states 
(2.40) to be eigenstates of the t^^\ij)^ G must be an eigenvector of t^'^\ii) with eigenvalue 
v^'^^jj). From (2.41) and (2.42) one can read off the eigenvalues of t^^-*(/u) 

r«M|A<", . . . , aW|G) = ^WM|A«, . . . , Ai')|G) = 

=1 a(//- A) ') ^.^1 a{\\ ' -fi)J 

The unwanted terms A^.^^ and A^^^ are computed in Appendix A and their cancellation 
(which ensures that the states (2.40) are eigenstates of the transfer matrix t^^\h) ) leads to 
the following set of Bethe Ansatz equations (BAE) for the first nesting 



=i«(a1^^-A,) fia(Ap)-A«) y kJ. 

i^k 



(2.45) 

The diagonahsation of r^'^^/d) can be performed by a third Bethe Ansatz. Prom the Yang- 
Baxter relation for r^'^\fx) one can derive the following intertwining relation 

r^'^ (A - /x) (T(f (A) ® if.)) = (7;(f (^) ® T(f (A)) r« (A - ;.) . (2.46) 

The components of (2.46) needed for the construction of an Algebraic Bethe Ansatz are 

D(2)(^)C(2)(A) = -^—^C^^){X)D('){f,) + ^[^~^\ c^^\f,)D('){X) 
ayjj, — \) a[\ — n) 

A^'\fi)C^'\X) = _J-_C(2)(A)A(2)(^) + ^[^~^l c^^\i^)A'~'\X) ^2.47) 
a[X — fj,) a[fj, — X) 

C(2)(A)C(2)(//) = C(2)(//)C(2)(A) . 
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As the reference state we pick 

|0)r=(5j , |0)(^) = ®^i,|0)i^) . (2.48) 

The action of the nested monodromy matrix Ti^^ (/j,) on the reference state |0)^^^ is determined 
by (2.43) and we find 

A(^\„)\0f^ = l[a{„-X^)\0f^ 



i^(^)(/.)io)(^) = n - A,) - bi^^ - A,)) io)(^) = n ^|x-i7yio)^'^ 



(2.49) 



We now make the following Ansatz for the cigcnstatcs of t*^^)(//,) 

|Af\...,A(f)) =C(^)(Af)) C(^)(Ar)...C(^)(A(f)) |0)(^) . (2.50) 

The analysis now proceeds analogous to the one for the first nesting. The eigenvalues of 
T^'^\lJi) on the states (2.50) are found to be 

r(^)M|Af',...,AW) = ^«M|Af',...,A(?)) = 



(2.51) 



under the condition that the spectral parameters are solutions to the BAE 

n''(^i"-^f)-n i' ■ p=i....,%. (2.52) 

i=l j=l «lAj, - Aj ) 

Using (2.51) in (2.45) {G is an eigenvector of t*^^)(A^"^'') with eigenvalue ^^''^•'(A^,^'') ) we obtain 
the BAE for the first nesting. If we then insert (2.51) in (2.44) and use the resulting 
expression for the eigenvalues of T^^\fi) in (2.34), we obtain the complete set of three BAE. 
If we shift the spectral parameters according to 

A. = A, + 1, aW^A^.^z, A(^) = A(^) + ^ , (2.53) 
these equations read 

l^k 

^e+^i I _ 3,(1) , _ 
jfc=l '^k 2 

Ap) -)^+^ 



Afc - 


A«- 


i 
2 


x(2) 


-A<" 


^ 2 


Am 


-A<" 


i 
2 


\(2) 
/\m 


-A<" 


i 
2 



n A; - Am ^ ^ _ -l-T -'j 2 , ^ 

T(2) _ r(2) _ . 11 S:(2) _ T(l) , i ' l,...,iV^. 
A^ Am t j=i Am Aj -p 2 
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Here we have used that n — + Ni — + + — N^, 77,2 = -^j,, which follow from 

our identification of the four (D^ basis vectors with the four possible configurations per site 
(1.1). The eigenvalues of the transfer martix r(/i) are computed by using only the wanted 
terms in (2.30) and (2.29), and inserting (2.44) and (2.51) in the resulting expression 

n ^) ^'"(.) . ( n ^) 
-^^^(^) - n n - \) -^^H/^) + n 77^^ (2.55) 

i=i «(/^ - >H ) \ 3=1 J i=i «(\ - f^) 

/ ( - \(1 

w - - 1 n n - A? •) + n n ^ 



1=1 - /i) j=i a{n - A} j=i - ^) 

Note that the condition of cancellation of the residues of the poles in (2.55) implies the 
BAE (2.54). The states (2.25) are simultaneous eigenstates of and t{^), and the energy 
eigenvalues can be determined from the eigenvalues of the transfer matrix (2.55) by using 
the trace identity (2.16). We find 

N,+Ni N,+Ni (2.56) 

3=1 3 + 4 3=1 

where we have reparametrised A^- = \cot{^). The eigenvalues of the momentum operator 
(2.18) can be obtained by setting /i = in the expression (2.55) for the eigenvalues of the 
transfer matrix, and then taking the logarithm. We find 

p(A„...,AJ = -^ ^ log = . (2.57) 

3=1 \^j~2j 3=1 



2.5. Algebraic Bethe Ansatz for the BBFF grading 

The construction of the NABA for this choice of grading follows roughly the same lines 
as for the FFBB case. Due to the fact that we are "truncating" fermionic lines, there are 
additional complications on the first two steps of the NABA. The situation is analogous to 
the one for the BFF grading in the supersymmetric t-J model, which was discussed in [9]. 
The necessary steps are all explained in that paper, so that we will constrain ourselves to 
only a brief sketch of the derivations, and mainly just quote the results. 

The grading is chosen in such a way that and 64 are fermionic (representing spin 
down/spin up electrons, respectively) and and 63 are bosonic (doubly occupied/empty 



14 



e2 = 0. The 



site) . In terms of the Grassmann parities this means that — — 1 and 
reference state in the k^^ quantum space jO);;. and the vacuum |0) of the complete lattice of 
L sites are fermionic with all spins up, i.e. 



|0). 








|0) = ®ti|0). 



The choice of grading leads to the L-operator 





a(A) + b{\)ei' 


bWel' 


bWef 


bWet' 






a(A) + b{X)ef 


b{X)ef 


b{X)ef 






b{\)el^ 


a(A) - b{X)ef 


-b^et;' 


\ 


bWel' 


b{\)ef 


-b{\)ef 


a(A) - 6(A) 



,44 



(2.58) 



(2.59) 



We again partition the monodromy matrix like in (2.21), which implies the following expres- 
sion for the transfer matrix 



(2.60) 



r^(A)|o) 



|o) 



(2.61) 



The action of the monodromy matrix on the vacuum can be computed using (2.58) and 
(2.59) 

f{a{X)f X 

{a{X)f 
(a(A))^ 

Vc,(A) C,(A) C,{X) (^)y 

As before we identify C„ as creation operators with respect to our choice of reference state 
and make an Ansatz of the type (2.25) for the eigenstates of the transfer matrix (2.60). The 
relevant intertwining relations are 



AM Ce(A) = (-l)^<.S+^«+^^ W(/^- A)5 C^{X) AM 

+ (_l)(ea + l)(e6 + l) ~ ^) 

1 



a(/u — A) 



^6(/«)^ac(A) 



(2.62) 



a(/j, — A) a{X — jj) 

CaS\) CaS^2) = r^FB(Ai - X.^Z ^,^^2) C,, (AJ 



where 



(2.63) 



Here ^ppB ^bbf the permutation matrices for the gradings = €3 = 1 , £3 = 
and = €2 = , €3 = 1 respectively. The intertwining relations (2.62) imply the following 
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action of the diagonal elements of the monodromy matrix on the states |A^, . . . , A„|F) 



k=l ^«l-«n . ^ 



(2.64) 



...ai 



(2.65) 



n ^ n 

n n 

+ E (/^) n s (A,) io)F«"-»^ , 

k=l 3=1 

where 

r^{i^:X = (-1)" ^^H/^ - aj--;l«,(/. - A„_,)t::r. • • - ajk.x ^^^^^^ 

X (-i)^<=Er=i'K+^+Er=i'^ciK+i) 

Here all indices Cj and c are summed over. The L-operators L^^\\) are constructed by 
multiplying the R-matrix r^^^(A) by the permutation matrix n^g^ and are found to be 

/a(A)+6(A)e^i h{\)ef h{\)ef \ 

L«(A)= h{\)e'^ a{X) + b{X)ef b{X)ef . (2.67) 

V biX)el' b{X)ef a{X)-b{X)efJ 

The expression for r^^^fj.) is significantly different from the one in the FFBB case, but by 
using an appropriate definition for the graded tensor product r^^^ (//) can again be interpreted 
as the transfer matrix of an inhomogeneous spin model on a lattice of n sites. Our reference 
state |0) is now of fermionic nature and we have to define a new graded tensor product 
reflecting this fact 

{Fmti = Fa,G^a (-l)(^^+^)(^''+^^) . (2.68) 

Effectively the new graded tensor product switches even and odd Grassmann parities, i.e. 
e„ — > + 1- terms of this tensor product the transfer matrix t^^\ii) given by (2.66) can 
be rewritten as 

r^'\l^)tX = MTi'\f^)) = ^tr(LW(/. - XJ^L^Un - A„_,)® . . . ®lS^V " Ai)) • 

(2.69) 

In the expression after the second equality of (2.69) we have explicitly written the tensor 
product (8) between the quantum spaces over the sites of the inhomogeneous model (the 
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L-operators are of course again multiplied as matrices). As before - must be an eigen- 
vector of T^^\ii), if . . . A„|F) is to be an eigenstate of t(//). The unwanted terms are 
computed in Appendix B, and the condition of their cancellation 

- ^"-"^ = (2.70) 

leads to the conditions 

F'^"-«^ = (a(-A,))^(rW(A,)F)""-"^ , k^l,...,n . (2.71) 

In the solution of the first nesting there is an additional complication compared to 
the FFBB case. Because of our change of tensor product, the L-operators L^^\X) are not 
intertwined by the R-matrix r(/x) defined in (2.63), but by the R-matrix 

rif^C = bWa, ^cci + a{y)6,, 6,, (-l)^^+^c+^„e, _ ^2.72) 

The following intertwining relation holds 

f(A - „) TW(A)®r«(/^) = rW(/^)®rW(A) f(A - n) . (2.73) 

If we make the following choice of reference state for the first nesting 

\0)i^= j , |0)« = ®Li|0)i'^ , (2.74) 

one can show that the model defined by (2.73), (2.74) and the graded tensor product ® is 
isomorpic to a model of the permutation type with BBF grading (describing two species of 
bosons and one species of fermions) and the graded tensor product (8). Effectively the first 
nesting thus describes the inhomogeneous model obtained by just "truncating a fermionic 
line" from the original model. This holds generally for u{m\n) symmetric models'^^'. 

In the next step of the NABA we thus face the task of solving the inhomogeneous BBF 
model. As in the first step we have to truncate a fermionic line, which can be done in a 
way analogous to the procedure described above. The intertwining relations arc of the same 
form as (2.62), but the R-matrices are now the ones of the BB and FF models respectively. 
Taking only the wanted terms into account we obtain the eigenvalues of the transfer matrix 
T^^) (jj,) as given below. The condition of the cancellation of the unwanted terms now reads 

n 

' k-l,...,n^ , (2.75) 

3=1 

where we have made an Ansatz of the form (2.40) for the eigenstates of the transfer matrix 
T^^\li), and where t^'^\ijl) is the transfer matrix of the second nesting given by the following 
expression 

r'^KA±, = - Am)r.;'.;:<-(A. - A« ,):;:::,:;:: . . . Lf\, - A«)-, . 

(2.76) 
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Here L^^\n) — a{^)I^'^^ + ^(/x)]!^^. Note that the roman indices in (2.75) and (2.76) can 
take the values 1,2, both of which carry bosonic grading and thus do not generate any 
minus signs. The second nesting can be shown to be equivalent to the inhomogeneous spin 
1 Heisenberg XXX model, or BB model in our terminology. The eigenvalues of r*^^)(/i) and 
the BAE for the inhomogeneous BB model are easily determined, and we can put together 
the three steps of the NABA like in the case of the FFBB grading. After reparametrising 
according to 

k = \-\. A« = A«-^, A(^^ = A(^)-^ , (2.77) 
the BAE take the form 

TT ^2 ^ TT ^ ^ +2 (2.78) 

k=l ^fc ■^j 2 m=l 2 

Ni r(2) r(2) . . 1(2) _ T(l) _ i 

r(2) _ r(2) _ - 11 T(2) _ T(i) i ' m-i,...,iV;. 
The eigenvalues of the transfer matrix are given by 




^=X/.-A«) 




- n 77T^ I 



A^i iV6 A^i 




fe=^i a(/^ - Af ^) ^ a{\f - 11) 

Energy and momentum eigenvalues are given by 

E\\,...X)-L- J2 I^ = -2 E cosik^) + L-2{N, + N^) 

j=l ^3 + 4 3=1 

p(\,...,AJ = -z 5] log(^y)= 5^ fc,. + (iV^ + l)7r , 
j=i + 2 j=i 



(2.80) 



where we have defined = ^tan{-^). Note that after the reparametrisations (2.77) and 
(2.53) the expressions for the BAE for the BBFF and FFBB gradings coincide up to the 
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substitutions <-> and <-> iV^. However there is no such obvious relation between 
the eigenvalues of the transfer matrix. For the energy eigenvalues it can be easily shown by 
taking logarithmic derivatives of the eigenvalues i^(^), that the expressions for BBFF and 
FFBB are identical up to the transformation A^^ ^ Ni and N-^ ^ apart from an overall 
minus sign. This is clear as under a particle-hole transformation for spin-up transforms 
into —H^. There is no such relation for the eigenvalues of the momentum operator. 

The BAE (2.78) are particularly useful, because they reduce to the BAE for the Suther- 
land solution of the supersymmetric t-J model in the sector without localons (A^; = 0). The 
classification of ground states and low lying excited states is particularly easy in this case. In 
our forthcoming publication we will analyse the bound state and ground state structure of 
this set of BAE, and describe the model by a set of integral equations in the thermodynamic 
limit [^^l 

2.6. Algebraic Bethe Ansatz for the FBFB grading 

Here the first step of the NABA is completely analogous to the one for the FFBB 
grading. The R-matrix and monodromy matrix of the nesting can be identified with the 
ones for FBF grading of the inhomogeneous supersymmetric t-J model, which was solved in 
[9]. The eigenvalues of the transfer matrix are found to be 

""""^ n ;;^( jg <.(.-v) (^-(.)- n ^] (2.1) 



Using the reparametrisation 



A, = A,+ |, A(^^ = A(^) + ^ , (2.82) 
the following form of the BAE is derived 

+ 2 / j=l \ - +2 

fe=l ^fe "t" 2 m=l "^j "I" 2 



Ni+Nt ~(2) _ ^(1) _ i 



j=l ^ 2 
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2.7. Algebraic Bethe Ansatz for the BFBF grading 

In this choice of grading we first truncate a fermionic hne analogously to the first step of 
the BBFF case. In the first nesting we truncate a bosonic line, which is done like in the first 
step of the FFBB case. The resulting expression for the eigenvalues of the transfer matrix 
are 



^'■'c) - n [ n - A.) j . n 
"'"w-n-^ n "(f-AS') 1- n 



t=la{^^-^k^) m=i \ j=i a{Xf^-fi) 

After reparametrising the spectral parameters according to 



A. = \-^, A(^) = A(^)-^ , (2.85) 



the BAE take the form 



j=l +2 
^b+Ni r ,(1) i TV, T(2) ,(1) 



TT I = n ^!!^ ^ I . = 1 AT, + AT, (2.86) 

A;=l 2 m=l ^"^ "T 2 

iV^+iV, ~(2) _ ^(1) _ i 

j=l Am Aj ~r 2 

The BAE (2.86) can be mapped onto the BAE (2.83) by a boson-fermion interchange like 
in the FFBB/BBFF case. 

2.8. Algebraic Bethe Ansatz for the BFFB grading 

The results for the eigenvalues of the transfer matrix and BAE are 
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L 



Ni+N^ ~ (1) 



K + I) j=i K - A^. +1 

^i+^i JV,+7V, r .(1) i T(2) ,(1) i 

/=1 '^j - '^l k=l '^k ~ '^j +2 m=l ~ '^j +2 

Af;+iV^ S^(2) _ _ i 



m=l '^"i '^j "T 2 



(2i 



where the we have reparametrised 



A. = A, + ^, A(^)=A(^)-^ . (2.89) 



2.9. Algebraic Bethe Ansatz for the FBBF grading 

This last possible choice of grading leads to the following results 



L N^+N, ~ (1) 



where 



,Afe + 1/ Afc - A^- +2 

A^^A^ A,-A^. -2 TT Am - A, - 2 

11 Ul) . 11 T . i 11 \(2) 

i=l '^j * jfc=l '^fc -r 2 m=l 2 

A^i+W, T{2) _ , (1) _ i 

(2.91) 

A.^A,-!, A(^) = A(^) + i . (2.92) 



Again there exsists a boson-fermion interchange map between (2.88) and (2.91). 
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3. Lowest Weight properties of the Bethe Ansatz states 



In this section we prove that all Bethe Ansatz states with finite spectral parameters 
are lowest weight states of the global u{2\2) symmetry algebra. For the case of the FFBB 
grading the lowest-weight conditions read 

0= r]\X,,...,XjF) = S\X„...,XJF) 
= Q,\X,,...,XJF) = qI|Ai,...,AJF) (3.1) 
= Q_,\X„...,XJF) = Qli|Ai,...,AjF) . 

Lowest weight theorems hold for all choices of the grading (the grading determines which 
generators annihilate the BA states), but without loss of generality we will perform the com- 
putations only for the specific case of the FFBB grading. The lowest weight property is of 
great importance as it first of all demonstrates that the Bethe Ansatz does not provide a com- 
plete set of eigenstates of the hamiltonian and because it can be used together with the 
global symmetry structure to generate a complete set of eigenstates. Similar computations 
have been performed for other models in [16-19]. 

We start by noting that in the basis of the Hilbert space introduced at the beginning of 
section 2.4, all lowering generators of the u{2\2) algebra are represented in the form X]fc=i ^f"- 
By using the correspondence of the four basis vectors over one site of the lattice with the 
four possible configurations (1.1), we find 



L 



e 



42 
k 



k=l k=l k=l /q r)\ 

L L L V-^-^J 

Q-1 = ^ ) Qi — ~ ^2 ' ^ ^ ' 

fe=l k=l k=l 

where e^-' are quantum operators acting nontrivially only on the /c*^ site of the lattice with 
matrix representation 

{eiT'-S,^5^, . (3.3) 

Using the representation (3.2) and the expression (2.20) for the L-operator we are able 
to rewrite graded quantum commutators as matrix commutators 

'(LMr , = (-l)(^^+^^)^^ ((e^^- LMY' - {LM e^^y) 

L } quantum \ (34) 

= (-l)fe+^.)- {[e\LM]^^,rJ' ■ 

where e^^ is a matrix of the form 

{e'r'-5^a8,, ■ (3-5) 
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Note that on the r.h.s of (3.4) e^^ and L^{fj,) are multiphed as matrices. The graded quantum 
commutator of with (T^(yu))"^ can now be expressed as a matrix commutator by using 
(3.4) and (2.21) 

L J quantum 

{{e'' LMT^'^' - (Lkif^) e^T^'^'') x (3-6) 
Summing over all sites of the lattice we arrive at 

iZ^^ATMr] = (-l)(^^+^^)- ib'^.TML^^J^'- (3.7) 

-k—i J quantum 

Equation (3.7) enables us to compute the graded commutators of the M(2|2)-lowering oper- 
ators with the creation operators C„(A). We find 



h C7„(A)] 


= KzD{X) - 


Aa(A) 




= K2D{\) - 




[Q-i,C,(A)} 


= WX) - 


AaW 




= ^aiC^iX) 




01i,C„(A)} 


= ^a2C3(A) 






= ^alC2(A) 





(3.8) 



We also note that all lowering operators annihilate the vacuum 

Q = n\0) = Q\\0) = QU\0) = Q,\0) = Q_,\0) = S\0) . (3.9) 



Lowest weight property for Q^, and rj 



In this subsection we show that all Bethe Ansatz cigcnstatcs are are annihilated by the 
m(2|2) generators = Q_i, J2 = Qi, and J3 = rj. Using (3.9) and (3.8) we find that 

Ja 1^1, • • • , A„|F) = 



i=l ) 

n k—1 n 

T.U^aS\)[Ja,C^,i\)} n C„^.(A,)|0)F'^»--H-1)^"^- 



j=k+i (3.10) 

J2U^aMi)i^aaD{X,)-A^^^{X,)) J] C„/A,)|0) F"--^(-l)^«n ' 



k=l i=l 

n k—1 



k=l i=l 

n k—1 



j=k+l 



=En^'',(^^) n ^^.,.(A,)io)(fi,)^-'"=-"^'=+-^"(-i)^»S'-^«< 

j=k+l 



k=l 1=1 
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The last equality follows by inspection of the r.h.s. of the previous hne, as D{X^) and 
■^aa^i^) be moved past the string of C^'s by using (2.26). Due to the cancellation of 
the unwanted terms (which are explicitly computed in Appendix A) all of the quantities Qj. 
vanish as we will demonstrate below. 

From (3.10) it follows that is the coefficient of the term, where the creation operator 
with spectral parameter is missing. It can be obtained from the term 



i=2 



(3.11) 



by moving {Sg^^D{X^) — ^aa^(Aj)) past YYj=2^a (^j) ^^^y taking "wanted" terms into 
account, i.e., by only using the first contribution of the intertwining relations (2.26) when 
moving ^^^(A^) or D{Xi) past Cj.(A^). The resulting expression for fl^ is found to be 



n 



«(Ai - A,) 



po„...ai 



-n 



«(A, - \) 
X {-l)^i=2%^a(^a{Xi))^ F""-'^i . 



(3.12) 



If wc compare (3.12) with the condition of the cancellation of the unwanted terms (•) and 
(•), we find that fli vanishes because 

(AiF)«''2-^ + (A^ir)»&2-b„ ^ Q _ 

The expression for Q,^ can be obtained in a similar fashion : 
We use (2.26) to rewrite \Xi . . . A^) as 



fc-i 



j=k+l 



i=l 



i=l 



\ck-2dk-l 



-2°'k-2 



(3.13) 

Now the action of the lowering operators J„ on the Bethe Ansatz states can be written in 
the form 



|Ai,---,A„|F) - 



fc-i 



Ja , ^'.^(A.) n^'',(^^) n ^a,(A,) s{x^':^:x \ io) F""-"^ 

j=k+l 



(3.14) 



1=1 
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Going through the same steps as in (3.10) it follows that can be obtained in essentially 
the same way as Jl^ with the result 



bi--bk-iah+i-K 



n ^ n 

n a(X -X,) n ^'^jb, ^'(AjaV.V.at'" 
1 ^ » j=k+l 



pa„...ai 



En 



X r{\k - r{K - KtlZ ■ ■ ■ r{K ' K-itllZ'-l 

x(a(A,))^ . 

(3.15) 

Inspection of the (•) and (•) in Appendix A shows that the condition of the cancellation of 
the unwanted terms 

implies the vanishing of flf.. This completes the proof of the lowest weight property of the 
Bethe Ansatz states for the generators Qi, Q^i and rj. 



Lowest weight property for Q\ and QLi 

In this subsection we show that the lowering generators Q^i annihilate all Bethe Ansatz 
states with finite spectral parameters A^ . We first reexpress in terms of matrix elements 
of the monodromy matrix (2.21) 

L 

Q-1 = ef = 1™ 

(3.16) 

< ' n—foo 
k=l 

Using the intertwining relation (2.26) we are able to determine the action of ^ ^a&(A*) 
states |Ai . . . A„) in the limit ^ —>■ oo 

n b b 

-i lim^A,,(A.)|A,...AjF) = -^nC,,(A,)|0) lim ^ (r«(^)'^y - " .(3.17) 

1=1 " 

This implies that 

Q\ \X,...XjF) = f[C,^{\M jim^/x(s«(^)F) 



1=1 



QU |Ai...AjF) = nC^6,(AJ|0) jim -^/.(5(^)(^)F) 



(3.18) 
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We now define new operators Ja^^ according to 



JW= lim -if,Bi'\f,) , a = 1,2 . 



ix-^oo 



(3.19) 



From (3.18) it is clear that tlie lowest weight property of the Bethe Ansatz states with respect 
to Qii is equivalent to the conditions 



J«F = 0, a = 1,2 



(3.20) 



RecaUing that 



"1 



F = J]cg)(A5^))|0)WG'S-^i 



we realize that we now face a problem very similar to the one we just solved for the case of 
T] and Qii. The graded commutator of j'^^ with C^\\'^^^) can be computed by using (2.37) 
and (3.19) 

'j':\c'^\>^''')\-8a,D^'\>^''')-A^^X^'^) . (3.21) 
The action of J^^ on the states jA^^'' . . . Xnl \G) now takes a form very similar to (3.10) 
J'^' |A?\...,A(;)|G) = 

fc-1 n-l 

T.Y[c'a^>^')[J':\c^:^>^t')] n c^a]i>^f)W' ^-.-(-i)-^^^^ 



fc=i j=i 

n, fc-1 



j=fe+l 



,.ai 



fe=l i=l 

X (-1)^"^ 

ni fc-1 



j=fe+l 



1 



fc=l i=l 



j=*;+l 



(3.22) 

The computation of the quantities Jl^.^^ follows the same strategy employed above for the 
determination of Q^ , and after some computations we find that all Q^^^ vanish, because the 
unwanted terms for the nesting cancel, i.e., 

This completes the proof of the lowest weight property for QLi- 



26 



Lowest weight property for S 
Noting that 

L 

S^y^ef^i \lvc.^xAM (3-23) 

we find 

-5 |Ai . . . AJF) = i lim AM\K ■ ■ ■ K\F) 

/i— >oo 

n b b 

' ^ ^ai-a„ (3.24) 

= inC,(A,)|0) ncl;'(AW)|0>C) 1™^(t<?M'-^g)°-'^'" . 

Z=l fe=l 

Therefore the condition that S annihilates all Bethe Ansatz states with finite spectral pa- 
rameters is equivalent to 

= z hm /.5(^)(^)|Af) . . . Ag) =: J^'^\X? . . . Ag) . (3.25) 
From the intertwining relation (2.47) we deduce 

[J(^),C(^)(A(^))] =D(2)(A(2))-^(2)(A(2)) , (3.26) 
which enables us to show by direct computation that 

j(^)|Af)...Ag)=0 
due to the cancellation of the unwanted terms for the second nesting 

(Af) + (A?^) = . 

This completes our proof of (3.1). 

4. Higher Conservation Laws 

In this section we derive explicit expressions for the conservation laws i/^g) and if(4) 
(which involve interactions between 3 and 4 neighbouring sites respectively). We use a 
generalisation of Tetcl'man's methodt^^'^^' to the supersymmetric case'^'. 

Let us define the "boost" -operator 
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where density of the hamiltonian given by the right hand side of (2.15). This 

operator obviously violates periodicity on the finite chain, but if used in commutators (which 
"differentiate" the linear n-dependence) it yields formally periodic expressions. 

The integrals of motion can be successively obtained by commutation with the boost- 
operator'^J 

where 

B^-Y,nHl^^, . (4.3) 

n 

Using (4.2) we obtain explicit expressions for higher conservation laws as follows : We 
have shown in [1] that the hamiltonian can be written in terms of -^(21 2) generators 
{ J,,J« = 1 ... 16} = {S^, S], SI ry„ ^J, ryl, J,, X„ Q^. g],,, Q^. g] ^ ±i| ^s 

i=i i=i i=i 

The structure functions of «(2|2) are defined by 

if(3) can be obtained by commutation with the boost operator B 

H^,)=t[B,H^,-^]=z[B,H'] 

L 

y^[-^fc+l,fc+2) Hk,k+l\ 



L 

I 

k=l 



(4.6) 



— i K"^ K^^ ffi^''Jk-l,a '^k,e '4+1,5 

k=l 

The next highest conservation law can be computed along similar lines and we find 



— ~2 K K K fp^ f^^ -4-1,0 4,e 4+1,0; 4+2,1/ 
k=\ 

L L 
_|_ pfe-l,A;+l _ 2 j];fc,fe+l 



(4.7) 



fe=i fe=i 

where p'^-i.'^+i is a graded permuation operator between the sites k — 1 and k -\- 1 with 
definition 

pk—l,k+l _ -Qfc— l,fcj-|-A;,fc+lj-|-fc— 1,A; 
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A. Appendix : Computation of the "unwanted terms" for the FFBB grading 



In this section we compute the so-called "unwanted terms" in the expressions (2.29) and 
(2.30). These unwanted terms are characterized by containing a creation operator with 
spectral parameter (SP) jj, in place of a creation operator with SP \^ (or A^^-*, A^-^-* for the 
nestings). The condition of cancellation of the unwanted terms leads to the Bethe equations. 
In order to obtain the expression for Aj., we first move the creation operator with SP Aj^. to 
the first place in (2.25), using (2.26) 

n k—1 n 

1=1 i=l j=k+l 

To get an unwanted term, we now have to use the second term in (2.26) to move D past 
Cfj^{Xj^), and then always the first term in (2.26) to move D (which now carries SP A^i.) to 
the very right, until it hits the vacuum, on which it acts according to (2.24). This way we 
obtain 



The computation of is more complicated. We first derive an expression for the con- 
tribution of A^^{iJ,), which we denote by A^_^. Proceeding along the same lines as in the 
computation of A^., we find 



1 (^iK - \) 

i^k 

X r{X, - X^t^ r(A, - X^f^ . . . r(A, - X,_^1:^ (A.3) 
X r(Afe - Afe+i)^^'^"J^'_^ r(Afc - Afe+sjb^+j^^' • • • - \)6^d„'_2 

x{a{X,)r5,^_^^A-ir^^'^''"' ■ 
The 6^ stems from the action of A^^^ _^ (A^) (which is what we get after moving A past 
all the C's) on the vacuum. We also had to include a 5^ ^ due to the fact that in (2.30) we 
denoted by the index of the C with SP fj,. The contributions of ^22 ( A*) ^33 ( A*) can be 
obtained in a similar way, and A^ = A^^^ -|- A^. 2 + Aj^. 3 is found to be 

X r{X,-X,)ti:ir{X,-XX2 ■ ■ A\ - K-XtlZZl i-^^^''^'"' 

X r{Ak '^k+l)bk+^d^-i ^y^k ^k+2)bk+2d^: ■■■^K^k ^n)b„d„_2 ■ 

(A.4) 
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This expression can be simplified by carrying out the contractions over the summation indices 
Ci, . . . ,c^. Noting that 

we are able to perform all c-summations with the result 

^(A.):;a <\ - • • • ^(^^ - ^^-1) t:2:: = n ^^^^ ^^.-.a, ■ (a.6) 

1=1 

Now we transform the remaining r-matrices into L-operators, using the identity 

rWf, = (r(A)nW)^,^ (-1)^"^= = {-iy^^^^L^'\XrJ . (A.7) 
Thus we obtain our final form for the unwanted terms due to —A^^lfj,) — ^22 (a*) + ^33 (a*) 

(A 8) 

We now insert (A. 8) and (A. 2) into the condition (2.33) for the cancellation of the unwanted 
terms, and multiply the resulting equation by the inverse of S{Xi^)a\',','\, which satisfies 
iS~^{Xi^))l^^"'l^ 5'(Afc)aj'.'.'.afe = Y['i=i^ai,Pi^ which is computed via (A. 5). After carrying 
out these manipulations wc arrive at (2.34). 

The unwanted terms for the nestings can be computed along similar lines, we refer to 
[9] for more detailed derivations and just give the results. For the first nesting we find 

^ ^ «(/^-Al^^)fia(AW-Af))}i 



and 



(A.9) 



(aWg)'^-'"^ =5«(A«)ft F^r.,.-9k,.fk-A - ^(\ -^M ]T J- , (A.10) 



where 

^^'HaI'V/;:/: = 



(A.ll) 
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The unwanted terms for the second nesting are 



,(2) _ 0[^^-K ) TT 1 TT„m(2) _ x(i)> 

^' -a(.-Ar)Ma(Af-Af)M ' ' ' ' 



(A.12) 



X(2) KAf-zO f^ 1 .^(Af - Aj") 



In all the above computations we have used that no two spectral parameters coincide, even 
if they belong to different nestings (e.g. ^ \/j,k). For the Bethe Ansatz without 
nesting the wave function vanishes if two spectral parameters coincide. This is sometimes 
called "Pauli principle in rapidity space" . The implications of coinciding spectral parameters 
in the nested Bethe Ansatz is currently under investigation'^^1 . 

B. Appendix : Computation of the "unwanted terms" for the BBFF grading 

In this section we compute the unwanted terms in the expressions (2.64) and (2.65). 
To obtain the expression for Af^F we first move the creation operator with SP to the 
first place in (2.25), using (2.62) 

n k—1 n 

''^FFsi^k ~ '^k-l)<:'k-ia-k-i ^ FFb{\ ~ \-2)ckJ>'^k-2 ■ ■ ■ F F b{\ ~ ^lYe^a^ 

We then use the second term in (2.62) to move D past (^^^(A^,), and then always the first 
term in (2.62) to move D (which now carries SP A^) to the very right, until it hits the 
vacuum, on which it acts according to (2.61). We find 

fA.FV"-'" = ^(A,)''-^^ ( KVlifl^ A 1 f <K) \ .32) 

V ' ^ '^"-'^'^ \a{K-lJi))\_\a{X,-X,)\a{-X,)) ^ ^ 

The computation of the unwanted term A^,F proceeds along the same lines as in the FFBB 
case (Appendix A). After some elemetary algebra we arrive at 

(A.^)'-'» - s(x^;:x (^<^) n ^^^^ M 

iT^k 

X r{X, - X,)ti:i r(A, - A^)?^ • • • r{X, - X,_,)tlr\Zll (B-3) 

X ^(Afc - Xk+i)bl^\dl-i ^'^^k - ^k+2)bllldl'^^ ■ ■ ■ ^iXk - K)bldl_2 

X (_i)'^fc+^^=i' ^;fc^'''*+^^(_l)(%+i)(%+i) 
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The summations over , . . . , e^^. can be carried out by using the identity 

Evaluation of the condition of cancellation of the unwanted terms (2.70) then leads to (2.71). 

C. Appendix : Equivalence of the BAE 

In this appendix we estabhsh the equivalence of the six different sets of BAE by means 
of particle- hole transformations in rapidity space^^l^. Let us denote the various sets of BAE 
by their underlying grading, i.e. FBBF, BFBF, BBFF, BFFB, FBFB, and FFBB. We first 
will demonstrate the equivalence of FBBF and BBFF. 

{i) FBBF ^ BFBF 

We start by expressing the third set of (2.91) 

^l+^l x(2) _ X (1) _ i 
j=l '^j "T 2 

as a polynomial equation of degree A^^ + A^ 

p(w) = n - - ^) - n - + = • (c-2) 

j=i j=i 

Among the A^ + A^j^ roots of (C.2) we consider A^^ roots w^, . . . , Wj^^, which we identify with 
\ . . . , The A; other roots of (C.2) we denote by w'j. Using the residue theorem we 
can derive the following equality 

where Cj is a small contour around Wj. The branch cut of the logarithm extends from 
— + f to 2^ = — |. By deforming the contours on the r.h.s. of (C.3) we arrive 
the the following equality 

g-'-(l^ )'-S-'-(l^)--(gS) ■ 

^We thank P.A. Bares for explaining this method to us. 
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where the last term on the r.h.s. comes from integration around the branch cut. The form 
of the polynomial p now implies that 



m=l 

N,+N^ 



(C.5) 



m=l 



Inserting (C.5) into (C.4) and exponentiating the result we obtain the identity 



n ^^ "^j ^ 2 _ TT fc 2 TT — Am + « 

>(1) i - 11 >(1) , i 11 >(1) >(1) 

j=l -^j - 2 k=\'^l ~ ^ fe + 2 m=l ~ 



Now we recall that Wj = Xj and use (C.6) in the second set of BAE in (2.91) with the result 

fc=l "^l ~^ 2 rn=l m i" 2 

This is precisely the second set of BAE for the BFBF grading (2.86), if we make the identifi- 
cation w\ = Xf\ The third set of the BAE (2.86) is also fulfilled by the spectral parameters 
Xf\ because they are roots of the polynomial equation (C.2). The first set of BAE for the 
BFBF and the FBBF gradings are identical anyway, so that we thus we have established the 
equivalence of the BAE (2.91) and (2.86). 



(ii) BFBF ^ BFFB 

We start with the first set of the BAE (2.86) for the BFBF grading 

^i+^i I >(i) 

,^fc^2/ 3=1 +2 

and rewrite it as the polynomial equation 



Partitioning the roots of (C.9) in the two sets {w^\k ^ 1 . . . N^+N^} = {X^\k ^ 1 . . . N^+N^} 
and {Wj\j — 1 . . .L — Nf^} and proceeding along the same lines as in step (i), we can derive 
the following equality 
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Recalling that — \ and inserting (C.IO) into the second set of BAE for the BFBF (2.86) 
grading and using that L — N^, + N^, we obtain 

r(2) x(l) i N,+N,7(l) , i Ni+N^^i) Wl) • 

y) .w,: = n fr-J-T n >a) ,a, . j - i....,a^+aV (cn) 

m=l '^"i ~ '^j +2 1=1 '^j ~ + 2 fc=l '^fe ~ ^3 ~^ 

Identifying w'j with A^- this imphes the equivalence of the BAE (2.86) and (2.88). 

{Hi) BFFB ^ FBFB 

The equivalence of (2.88) and (2.83) can be shown by an analogous computation to the 
one in step (i). 

(iv) BFBF ^ BBFF 

We start by converting the second set of the BFBF BAE (2.86) into the polynomial 
equation 

k=l j=l k=l 3=1 

= . 

(C.12) 

We partition the roots of (C.12) into the sets {w,^\k = 1 . . . A^^+A^,} = {A^.^^ | A; = 1 . . . AT^+A^,} 
and {w'j\j — 1 . . .N^}. Proceeding along the same lines as in step (i), we can derive the 
following equalities 

iC 13) 

j—]^ Am ~ "^j ~r 2 1=1 I ~ ~ j=l ~ j ~ 2 

Inserting (C.13) into the BFFB BAE (2.88) and identifying the w'^ with Xf^ in the BBFF 
BAE, we see that (2.88) and (2.78) are equivalent. 

(v) FBFB ^ FFBB 

This step is analogous to step (iv). 
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